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ABSTRACT 

Galaxies above redshift 1 can be very clumpy, with irregular morphologies dominated by star complexes 
as large as 2 kpc and as massive as a few x 10 8 or 10 9 M Q . Their co-moving densities and rapid evolution 
suggest that most present-day spirals could have formed through a clumpy phase. The clumps may form 
by gravitational instabilities in gas-rich turbulent disks; they do not appear to be separate galaxies 
merging together. We show here that the formation of the observed clumps requires initial disks of gas 
and stars with almost no stabilizing bulge or stellar halo. This cannot be achieved in models where disk 
galaxies grow by mergers. Mergers tend to make stellar spheroids even when the gas fraction is high, 
and then the disk is too stable to make giant clumps. The morphology of high-redshift galaxies thus 
suggests that inner disks assemble mostly by smooth gas accretion, either from cosmological flows or 
from the outer disk during a grazing interaction. 

Subject headings: galaxies: formation — galaxies: high-redshift — instabilities 



1. INTRODUCTION 



In the standard A-CDM cosmology ( Blu menthal et all 
1984 J), galaxies as semble through hierarchical merging 



(Da vis et al.1 Il985l ) and then evolve under a combina- 
tion of internal and environmental processes. Most of a 
galaxy's mass comes from mergers with m uch smaller ob- 
jects and major merger s are less frequent (|De Lucia et al.l 
12001 iGenel etalEom . Major and minor mer gers tend to 
transform disks into sphero ids and ellipticals (|Naab et al.l 
l2007t iBournaud et al.ll2007[) . but if the gas fraction is high 
enough, then massive and extended disk structures can 
persist (|Robertson et al.ll2006t ISpringel fc Hernqu ist 2005; 
iRobertson fc Bullock! 12008V Therefore it may be possi- 
ble, within the hierarchical framework, to explain the for- 
mation of disk galaxies with morphologies and kinem at- 
ics like those observed locally (|Governato et al.1 [2007). It 
has recently been suggested, however, that a large part of 
the galaxy mass could instead come f rom diffuse gas ac- 
cretion, in particu l ar along cold flow s (|Dekel et al. 2009a; 
lOcvirk et ai]|2008t iKeres et al.ll2008th Thus, the build-up 
of high-redshift galaxies may result from two processes, 
one hierarchical and the other somewhat smooth. This 
Letter considers the resolved properties of these galaxies 
to suggest that the smooth process dominates in a high 
fraction of cases. 

Galaxies are inc reasingly clumpy with redshift 
(|Conselice et alj [2005h These clumps are not just features 
in ot h erwise normal spiral s and ellipticals ( Cowie et al.1 
119961: Ivan den Bergh et al.l Il996h . Most high-redshift 
galaxies do not have spirals in restframe blue and uv 
bands, nor do they have bulges or exponential profiles. 
Usually, a large fraction of their optical light (up to 50%) 
and luminous m ass (up to 30%) is confined to a few 
kpc-size clumps (Elm egreen fc Elm egreen 20051 hereafter 



EE05). Even the inter-clump l ight does not follow a 
spiral- or exponential- like profile (|Elmegreen et al.l l2~005). 
Highly aligned clu mps have been called chain galaxies 
(|Cowie et al.1 Il996l) . while rounder system s have been 
called clump- clusters (|Elmegreen et al.ll2004l) . Both types 
could be progenitors of modern spiral disks, viewed with 
different orientations (Elmegreen et al. 2005; Bournaud, 
Elmegreen & Elmegreen 2007, hereafter BEE07). If the 
clumps are massive e nough, then th e y spiral to the cen- 
ter to make a bulge (jNoguchil 119991: llmmeli et all 12004 
Elmegreen, Bournaud & Elmegreen 2008). 

Clumpy galaxies are so frequent at 1 < z < 5, and 
the clumps evolve so quickly, that most present day spi- 
rals could have had a clumpy phase in their past. Their 
comoving space density of the UDF is ~ 10~ 3 Mpc~ 3 be- 
tween z — 1 and 4, comparable to the density of spirals in 
the same r edshift range and with the same absolute zsso 
magnitude (jElmegreen et al.ll2007l ) . Moreover, the lifetime 
of the clumpy phase is short, ~ 0.5 Gyr, according to sim- 
ulations (BEE07) . If wc consider the ratio of the lifetime 
of the clumpy phase to the Hubble time as a measure of 
the fraction of galaxies in that phase at any one time, and 
multiply the inverse of this ratio by the observed space 
density of the clumpy systems, then we get a total space 
density for galaxies that ever went through the clumpy 
phase. This total is comparable to the space density of 
modern spirals. 

Several observations suggest that clumps in most chain 
and clump-cluster galaxies formed inside their disks rather 
than entered from outside in a merger. First, the distri- 
bution of the ratio of axes for the combined population is 
approximately flat, suggesting that most of the clumpy 
types are disks viewed in random orientations (EE05). 
Second, the clumps in chains are highly confined to the 
average midplanes, which makes external capture unlikely 
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(Elmegre en fe Elmegreenl l2006al ) . Third, the masses and 
sizes of the clumps are consistent with their formation by 
gravi tational instabilities for the observed velocity di sper- 
sion ([Forster Schreiber et a l. 2006; Genz el et al.ll2008l ) and 
a gas column densi ty comparable to the t otal in today's 
inner spiral disks ([Elmegreen et alj |2009[ ) . Fourth, the 
height of a clump is comparable to the disk half-thickness, 
suggesting th at both are determined by the gravitational 
scale length ([Elmegreen fe Elmegreenl l2006a|) . Fifth, the 
distribution of relative clump position in a UDF chain is 
the same as for edge-on clump clusters (Elmegreen 2009). 
Sixth, the largest clumps have similar ma sses and ages, un- 
like t he expectation for random capture ([Elmegreen et all 
2009). The shape of clumps is another clue to their in situ 
origin: they are not elongated like spiral arms, so they 
have to form quickly in a highly unstable disk (Toomre 
Q < 1). Red interclump colors and a monotonic rise in 
the r otation curve for one studied case (|Bournaud et alj 
2008) also suggest the clumps are part of a disk. 

Here we show that giant star-forming clumps require 
high turbulent speeds and a dense disk with few stars in 
a spheroid. Simulations of galaxy mergers predict some- 
thing different: that a large fraction of the stars should 
end up in a spheroid. Such a spheroid stabilizes the disk 
and prevents the observed clumps from forming. 

2. REQUIREMENTS FOR GIANT CLUMP INSTABILITIES IN 
PRIMORDIAL GAS DISKS 

2.1. High velocity dispersions 

Typical attributes of UDF clumpy galaxies are a stellar 
mass M* = 6 x 10 10 M and a disk radius R — 9 kpc 
(EE05). The gas mass fraction / g (gas to total baryon ra- 
tio) is not known o bservationall y. An estimate of / g ~ 
50% was made by IDaddi et all (|2008l ) in BzK-selected 
galaxies. We apply this fraction to clumpy galaxies be- 
cause (i) Daddi et al.'s galaxies are clumpy in the HST 
images, and (ii) other BzK-selec ted galaxies have mas- 
sive clumps like the UDF g alaxies dForster Schreiber et alj 
I2006t iGenzel et al.ll2008j ). iBouche et al.1 (|2007f ) also sup- 
port gas fractions around 50%. Much larger gas frac- 
tions w ould exceed the dynamical mass from cir cular ve- 
locities (|Bournaud et al.l l2008t IDaddi et al J 12008( 1. Lower 
gas fractions would make the clumps too big for the ob- 
served velocity dispersion. Thus, we assume f s = 50% 
inside the stellar disk radius, and a gas mass M g = M* = 
6 x 10 10 M Q . 

The Jeans length for gravitational instabilities is Aj = 
a 1 1 (7tGS). For somewhat uniform disks, £ ~ M g / (ttR 2 ), 
which gives, 

We consider clumps of size Aj ~ 500 pc or larger, so 
a ~ 50 km s~ 1 for the gas that forms giant clumps. Lower 
velocity dispersions would give smaller clumps. This result 
is consistent w ith Ha observations that sugg e st a high tur- 
bulen t speed ( Forster Schreiber et al.ll2006t IGenzel et al.l 
200 8j), and with observations o f moderately thick disks 
(|Elmegreen fc Elmegreenll2006al l. 



2.2. Dense disks with low-mass stellar spheroids 



2.2.1. Analytical constraints on Q 

Clumps form by gravitational instabilities if Q = 
an/ (ttGY.) < 1. We consider first this requirement on 
Q for a model in which (1) M g = 6x 10 10 M in a disk of 
radius R = 9 kpc with a turbulent speed a — 50 km s _1 ; 
(2) the total stellar mass of M* = 6 x 10 10 M Q inside R 
is partly in a disk of mass M*^, partly in a central bulge 
of mass M*^ and radius 1 kpc, and partly in a stellar 
spheroid with M*,jj=1-M*,£)-M* i _b; (3) there is always a 
dark halo with a mass /d x (Mg+M*) inside R. We choose 
/d=0.5, which means that 2/3 of the mass inside R is bary- 
onic and 1/3 is dark. High redshift kinematic observations 
(Daddi et al. 2008, Bournaud et al. 2008) are consistent 
with this fraction. Modern disk galaxies should be about 
the same because late infall of dark matter and baryons 
typically follow each other (Semelin & Combes 2005). Lo- 
cally, the dark-to-disk ratio fr> is on average 0.6-0.7, and 
most generally in the 0.3-1.0 range (e.g., Persic & Salluci 
1990); our frj — 0.5 choice is thus realistic. We further 
discuss the influence of this parameter in Section 3. 

We compute Q for the combined gas and stellar disks, 
assuming the stars have about the same velocity disper- 
sion as the gas because both are heated by gravitational 
insta bilities and clump interactions. Then (W ang fc Silkl 
11994 : 

Q(r)~f^^ + 3 - 36G5 )«f r V 1 . (2) 
\a gas n(r) o- Btar K{r) J 

When a is the same for gas and stars, Q ~ 
<jk/ (7rG [£ gas + £ star ]). Profiles of Q{r) are shown on 
Figure 1. When all the stars are in the disk, £ is high 
enough and n is low enough that Q < 1 over a large part 
of the disk. When less than 80% of the stars are in the disk 
and more than 20% are in a halo or bulge, Q > 1 over the 
whole disk. Then spiral arms can form with star-forming 
clumps inside, but isolated and round clumps become less 
likely. 

2.2.2. Numerical Simulations 

To check these simple Q estimates, we ran models of 
galaxies with various ratios of disk: bulge: spheroid stellar 
mass and determined when giant clumps appeared. The 
mass and size parameters are as in the calculations above, 
and we start with cr=50 km s _1 for gas and stars. The sim- 
ulations were run with higher resolution than in BEE07: 
3 x 10 6 particles each for the stars, gas, and dark matter, 
and a spatial resolution of 30 pc over the whole disk. The 
Jeans length is more than ten times larger than the reso- 
lution in our initial set-up, so numerical fragmentation is 
avoided (jTruelove et al.lll997[ l. 

Figure 2 shows the mass distribution of gas and stars 
at the most clumpy instant for several models. When 
100% percent of the stellar mass is initially in the disk, 
dense round kpc-sized clumps form as in the observed high- 
redshift galaxies. When 80% of the stars are in the disk 
and 20% are in a bulge or extended stellar halo, the disk 
is still somewhat clumpy, but the clumps are not isolated, 
round and gravitationally bound - they are bright spots 
in shearing spiral arms. When 60% or less of the stars lie 
in the disk and the rest are in a bulge or halo, the disk 
forms mostly spiral arms (small, low-mass clumps can still 
form, but only when a is small). 
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These models confirm the earlier calculations based on 
Q alone. The formation of giant clumps requires a and £ 
to be high and k to be low (Eq. 2). Stars in the rotating 
disk contribute to £ and help drive the clump-forming in- 
stability, while stars in the bulge or spheroid increase k and 
stabilize the disk. Thus the formation of giant clumps in 
galaxies with masses and sizes typical of the UDF clump 
clusters requires that almost all of the stars lie in a ro- 
tating disk together with the gas when clump formation 
starts. If this constraint is not satisfied, then Q > 1 and 
massive clumps do not readily form. If turbulent dissipa- 
tion decreases a enough to form clumps when the stellar 
disk mass is relatively small, then the resulting clumps 
will have low masses, unlike the observed clumps. Such 
small clumps are seen in the outer parts of the 60:30:10 
and 60:10:30 models in Fig. 2. The large observed sizes of 
chimps at high redshift imply that a is high, and values of 
cr higher than the one conservatively assumed above would 
make the formation of clumps even harder to achieve if a 
stellar spheroid is present, by increasing the Q parameter. 

3. CLUMPY GALAXIES IN THE HIERARCHICAL CONTEXT 

The clumpy galaxies in the UDF are not obviously merg- 
ers, although outer tidal tails and other merger debris 
could be too faint to see. They are also not former mergers 
because most of their stars should be confined to a disk in 
order to get the giant clumps. An important characteristic 
of disks built by galaxy mergers is that they are embedded 
in stellar spheroids. For inst ance, the model spiral galaxy 
formed in a major merger bv lSpringel fc Hernquistj ([2005) 
has only ~ 25% of its stellar mass in the rotating disk, 
while ~ 50% is in a massive, kpc-sized bulge and ~ 25% 
is in a more extended diffuse stellar halo. Even mergers 
starting as pure gas rarely en d-up with more than 50 % of 
the stellar mass in the disk: iRobertson et alj |2006) find 
that for gas fractions of 0.4 or 0.6, equal-mass mergers 
end up with only 10% or less of the stars in the disk, the 
rest being in a bulge and a stellar halo; even initial gas 
fractions of 80% produce disks containin g only 40-45% of 
the s tellar mass (models DC, EC, FC in IRobertson et ail 
2006). In the case of minor mergers, their models DCm, 
ECm, FCm indicate final disk fractions of 60-75% after a 
single 8:1 merger, for gas fractions rang ing from 0.4 to 0.8. 
A larger study by ([Hopkins et al.ll2008f ) also finds that the 
fraction of stars in the disk after a single 8:1 minor merger 
is 70% to 85%, depending on the gas fraction and other 
parameters (e.g., their Fig. 12); after a 2:1 merger, the 
mass in spheroids is equal to or larger than the stellar disk 
mass. A 1:1 merger ends-up with less than one third of 
the stellar mass in the rotating disk component. 

If we consider a typical clumpy galaxy in the UDF and 
assume that, over the last doubling of its mass, mergers 
were the dominant growth process, then that growth may 
have occurred, for example, by a 1:1 major merger or by 
six successive 8:1 minor mergers. For gas frac tions of 50% 
in th e progenitor galaxies (likely for z ~ 2, iDaddi et alj 
2008j) , such mergers will leave only 20% or less of the stel- 
lar mass in the rotating disk, the vast majority being in 
the bulge and halo. Even higher gas fractions will not al- 
low more than ~ 50% of the stellar mass in the disk. If we 
assume that the last mass doubling event was half done 
by mergers and half by smooth accretion of cold gas, then 



this requires one 2:1 merger, or 3 to 4 mergers of mass 
ratio 8:1. These will leave only 40-50% of the stellar mass 
in the disk, given the various results mentioned above. If 
we consider mass assembly over longer periods, then the 
effect of mergers on the bulge and stellar halo will be even 
more dramatic. 

The ubiquity of giant clumps in z ~ 2 disk galaxies can 
be explained most easily if only a small fraction of the 
stars lie in a bulge or halo before the clumps formed. We 
found that this fraction should be < 20%, assuming a dark 
to baryonic ratio fr> of 50% inside the optical radius. This 
ratio ju should not be much lower in the progenitors of 
present day spirals. For a very low fu of 25% (80% of 
the mass inside R is baryonic), a slightly more massive 
stellar spheroid could be present without preventing the 
formation of large clumps. Still, the fraction of baryons 
that can be in a bulge or halo remains limited to < 30% 
even in this extreme case (to keep the total dark+stellar 
spheroid constant). This would allow a relatively higher 
contribution from mergers in the formation of these galax- 
ies, but it still implies that most of the mass has to enter 
in a smooth accretion process instead of through major 
and minor mergers. 

The key mechanism is unlikely to be distant interactions 
that de-stabilize the disk and provoke internal clump for- 
mation. Interactions can destabilize a disk (e.g., di Mat- 
teo et al. 2008) when the disk was marginally stable, 
but not when Q was originally high because of a mas- 
sive spheroid. Furthermore, an interaction-induced pro- 
cess would lead to centrally-concentrated star formation 
because of the associated angular momentum redistribu- 
tion. This would contradict observations often showing 
the star-forming clumps at large radii or in ring-like struc- 
tures in bulge-free disks (Elmegreen & Elmegreen 2006b; 
Genzel et al. 2008; Elmegreen et al. 2009). 

Gas physics and thermal cooling can affect the long-term 
evolution of clumps, as studied with various hydrodynamic 
codes by Immeli et al. (2004), Debasttista et al. (2006) 
and Tasker & Bryan (2008). However, in both our sim- 
ulations and our analytical calculations, thermal cooling 
is not a critical issue because the turbulent motions in 
the gas greatly exceed the thermal motions on kpc scales 
where the instability operates. For our simulations, the 
turbulent speed is set to a realistic value considering the 
observed clump sizes, disk thicknesses, and gas velocity 
dispersions. The medium in which we form clumps thus 
has realistic properties. Only the small-scale properties of 
the clumps, such as core formation, molecule formation, 
and dense cluster formation, all of which are unobserved 
so far, can be affected by thermal cooling. In a super- 
sonic medium, thermal cooling rates affect the thickness 
of the shock fronts, but not the overall dissipation rate 
of the turbulent motions. In this regard, we note that 
large-scale simulations of supersonically turbulent gas and 
clump formation require numerical methods that can han- 
dle high velocity dispersions in a thermally cold interstellar 
medium. 

If most of the mass assembly of these galaxies came from 
major and minor mergers, then the disk density would be 
too low and the shear rate would be too high. Clumps 
could not form with the masses and sizes typically ob- 
served at high redshift. For instance, the merger-produced 
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disk in the model by iRobertson fe Bullock! (2008) is mas- 
sive and extended, but it only grows spiral arms without 
giant clumps. The assembly of high-reshift galaxies, with 
frequent clumpy morphologies, is unlikely to be mostly 
driven by hierarchical merging of smaller galaxies. Even 
if minor and major mergers were responsible for only half 
of the growth, the spheroids would be too massive for disk 
clumps. The formation of giant clumps points to massive, 
highly turbulent disks that have relatively small bulges 
and stellar haloes. This requires that the dominant pro- 
cess of mass assembly be some smooth accretion of cold 
and diffuse gas. This conclusion is consistent with the re- 
cent picture in which young thick disks form by cold flows 
(Dekel et al. 2009a; K eres et al. 2008) and oth er types 
of diffuse gas accretion (jSemelin fe Combes! 2005), bulges 
form by internal, clump-d riven evolution ( Elmegreen et all 
120081: iGenzel et all 12008!) . and the thin disk forms later 
by fu rther smooth accretion (e.g., iBournaud fc Combes! 
2002). Our result does not imply that mergers do not 
occur at high rcdshift, but only that they cannot be the 
main mechanism for disk assembly. 

4. CONCLUSION 

The ubiquity of giant clumps in high-redshift disk galax- 
ies constrains the mass distribution and therefore assembly 
process. Most of the stellar and gas mass should be in a 
disk, rather than in a bulge or spheroidal halo, in order 
to get the clump masses, sizes, and morphologies correct. 



This implies that the assembly is mostly smooth, with 
only a small fraction through minor and major mergers. 
Bulges then form as an aftermath of clump evolution. An 
exception might occur for SO galaxies, which have much 
larger stellar spheroids; mergers could have played a more 
important role in their formation and structure. 

Cosmological models including warm dark matter could 
lead to the required smooth gas accretion because warm 
dark matter has relatively little substructure. The model 
by Heller, Shlosman & Athanassoula (2007) had a mono- 
lithic collapse of gas inside a single halo without hierar- 
chical merging: the gas flow was smooth, and the disk 
formed through a clumpy phase. Cold dark matter mod- 
els may also produce smooth gas accretion and satisfy 
our constraint on the disk mass fraction (e.g., Dekel et 
al. 2009a; Semelin fc Combes 2005; Keres e t al. 2 008). 
Recent models by iDe kel. Sari fc Ceverino I (|2009bh and 
lAgertz. Tevssier fc Moore I (|2009l) suggest that cold flows 
could maybe lead to the formation of clumpy galaxies. 

Present-day spiral galaxies have apparently evolved 
through a clumpy phase, and this observation may prove 
to be a key factor in understanding galaxy assembly and 
the nature of dark matter. 

We are grateful to Debra Meloy Elmegreen for useful 
discussions, and to an anonymous referee for constructive 
comments. Numerical simulations were carried out on the 
NEC-SX8R vector computer at CEA/CCRT. 
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Fig. 1. — Toomre Q parameter for the gaseous and stellar disks, with several assumptions about the stellar mass distribution. The legend 
indicates the proportion of disk:halo:bulge stars. We assume a uniform surface density for the disk, a Hcrnquist profile with a scale-length 
of 600 pc for the bulge, and the same profile with a 6 kpc scale-length for the stellar halo. There is always a dark matter halo, with a core 
Burkert profile (a cuspy profile would somewhat increase Q). 
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Fig. 2. — Surface density maps for the gas (top) and stars (bottom), at the most "clumpy" instant seen face-on. The stellar mass distribution 
in the disk:bulge:halo components is varied as in Fig. 1; there is always a gas disk and a dark matter halo. All snapshots are shown with the 
same logscale colorbar. 



